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We consider model A dynamics for a heating quench from the disordered (confined) into the ordered (deconfined) 
phase of SU(3) lattice gauge theory. For 4 N% lattices the exponential growth factors of low- lying structure function 
modes are calculated. The linear theory of spinodal decompositions is compared with the data from an effective 
model and the Debye screening mass is estimated from the critical mode. Further, the quench leads to competing 
vacuum domains, which make the equilibration of the QCD vacuum after the heating non-trivial. We investigate 
the influence of such domains on the gluonic energy density. 



1. Introduction 

In Ref. |1I2| it is argued that a heating quench 
of QCD from its confined phase of disordered 
Polyakov loops into its deconfined phase of or- 
dered Polyakov loops leads to vacuum domains 
of distinct Z3 triality, and one ought to be con- 
cerned about non-equilibrium effects due to heat- 
ing of the system. This comes because the heating 
is a rapid quench and the QCD high temperature 
vacuum carries structures which are similar to the 
low temperature phase of analogue spin models. 

Here we report a similar investigation for 
SU(3) lattice gauge theory and preliminary re- 
sults about the influence of such domains on the 
gluonic energy density and pressure of pure SU(3) 
lattice gauge theory. The Markov chain Monte 
Carlo (MC) process provides model A (Glauber) 
dynamics in the classification of Ref. 3 . As time 
step a sweep of systematic updating with the 
Cabibbo-Marinari 0] heat-bath algorithm and 
its improvements of Ref. |5I6| is used (no over- 
relaxation, to stay in the universality class of 
Glauber dynamics). Although this is certainly 
not the physical dynamics of QCD, in the present 
state of affairs it appears important to collect 



qualitative ideas about eventual dynamical ef- 
fects. For this purpose the investigation of any 
dynamics, which actually allows for its study 
ought to be useful. 

2. Preliminaries 

We report numerical results for the structure 
factors (or functions) defined by 



S(k,t) = (u(k,t)u(-k,t) / 
i(f,t)u(r\t) \ exp (ik ■ (f — r'^j drdr'. (1) 



Here u(r, t) is the relevant fluctuation about some 
average. For gauge systems we deal with fluc- 
tuations of the Polyakov loop (for analogue spin 
systems with fluctuations of the magnetization). 
The time t corresponds to the dynamical pro- 
cess, i.e., in our case the Markov chain MC time. 
Within the linear theory the differential equation 



du(k',t) 
dt 



w(fe') u[k',t) = g(k') 



(2) 



with g(k) = — r /g y exp(ifc • r) 
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can be derived, where T is the response coeffi- 
cient and /q the derivative of the coarse-grained 
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free energy density with respect to the fluctuation 
variable. The general solution of Eq. J2J is 



ii(k',t) = Cexp (u(k')t\ 



with = -r (Kk 2 + f 



(3) 
(4) 



Where X is the lowest order expansion coefficient 
of the free energy density about its minimum /o, 
i.e. K > 0. If /o > holds, Eq. J3J implies 
that the amplitude of any fluctuation approaches 
a constant exponentially fast with time. But if 
the second derivative is negative, then one sees 
an exponential growth of the fluctuations for mo- 
mentum modes smaller than the critical value 



\k\ < k c 



K 



1/2 



(5) 



The equation of motion for the structure factor 
Q is derived by taking the time derivative of 
S(k, t) and using Eq. 



dSjk, t) 
dt 



2u(k,t) u*(k,t)u>(k) 



+ ([u*(k,t) + u(k,t) g(k)) . (6) 



The average of fluctuations about the mean of the 
fluctuation variable has to be zero (u(k,t)^ = 0. 
Thus © becomes 



dS(k,t) 
dt 



2w(k){>(k,t). 



with the solution 
S(k, t) = S(k, t = 0) exp ( 2u(k)t 



(7) 



(8) 



Again, for f < low momentum modes grow 
exponentially. The value of the critical momen- 
tum is the same as for the fluctuations. Originally 
the linear theory was developed for model B |7I8| . 
Details for model A can be found in Ref. [2]. 

During our simulations the structure functions 
are averaged over rotationally equivalent mo- 
menta and the notation S ni is used to label struc- 
ture functions of momentum 

2-7T 

k = — n where \n\ — rij . (9) 
L 



We recorded the modes (including the permuta- 
tions) m: (1,0,0), n 2 : (1,1,0), n 3 : (1,1,1), n 4 : 
(2,0,0), n 5 : (2,l,0),n 6 : (2,l,l),n 7 : (2,2,0),n 8 : 

(2.2.1) and (3,0,0), n 9 : (3,1,0), n w : (3,1,1), 
n n : (2,2,2), n 12 : (3,2,0), n 13 : (3,2,1), n u : 

(3.2.2) , n 15 : (3,3,0), n w : (3,3,1), n l7 : (3,3,2), 
nig: (3,3,3). Note that there is an accidental 
degeneracy in length for n^. 

3. Numerical Results 
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Figure 1. The first structure function mode for 
pure SU(3) lattice gauge theory on 4 lattices. 

In Fig. the time evolution of the first struc- 
ture function mode after a heating quench from 
(3 = 5.5 — > 5.92 in pure SU(3) lattice gauge theory 
is depicted. For the 4 x 16 3 lattice the pseudo- 
transition value is (3 C = 5.6908 (2) |S] and for the 
larger lattices the value is expected to be close-by. 
Our results are averages over 10,000 repetitions 
of the quench for the 4 x 16 3 lattice and 4,000 
repetitions for the 4 x 32 3 lattice. The 4 x 64 3 
lattices are still running: Presently there are only 
60 repetitions and only some structure function 
maxima are presently reached. Notable in Fig. ^ 
is the strong increase of the maxima S™ x with 
lattice size. In our normalization non-critical be- 
havior corresponds to a fall-off ~ l/A" 3 and a 
second order phase transition to a slower fall-off 
- with < x < 3. As the N a -> oo limit 

is bounded by a constant, our figure shows that 
with our lattice sizes the asymptotic behavior has 
not yet been reached. 

To determine the critical mode fc r we fit the 
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Figure 2. Determination of k c for the pure SU(3) 
lattice gauge theory on 4 N% lattices. 

initial increase of the structure functions. The 
results are combined in Fig.|2and indicate ak c = 
2nn c /N<j « 0.34 (from the figure n 2 JNl sa 
0.003). Using results of Miller and Ogilvie [TU] . 
we have m£> = \^3k c , where tud is the Debye 
screening mass at the temperature in question. 
The relation k c /Tf — N T ak c , where Tf is the 
final temperature after the quench, allows us to 
convert to physical units. For our quench we have 
T f /T c = 1.57 and get 



m D = V3N T ak c T f = 3.7T C . 



(10) 



For pure SU(3) lattice gauge theory T c = 
265 (1) MeV holds, assuming a = 420 MeV for the 
string tension, while for QCD the cross-over tem- 
perature appears to be around T c w 165 MeV, see 
[I3j for a recent review. 

In Fig. ^ we observe that not only the height 
of the peaks increases with lattice size, but also 
the time i max , = S ni (t max ), which it takes 

to reach them. Whereas 5'™ ax has finally to ap- 
proach a constant value, £ m ax is expected to di- 
verge with lattice size due to the competition of 
vacuum domain of distinct Z3 triality. 

Using Fortuin-Kasteleyn (FK) clusters, the 
competition of such distinct vacuum domains can 
be made visible for the analogue Potts models [2J . 
The states of the three-dimensional, three-state 
Potts model substitute then for the Z3 trialities 
of SU(3) lattice gauge theory. In Fig. we com- 
pare the evolution of geometrical and FK clusters 
for a quench of this model from its disordered 
into its ordered phase. We plot the evolution of 
the largest clusters for the three Potts magnetiza- 
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Figure 3. Largest geometrical and FK clusters 
for the 3D 3-state Potts model quenched from 
= 0.2 to Pf — 0.3 at zero field on a 40 3 lattice. 

tions in zero external magnetic field h. While the 
the system grows competing FK clusters of each 
magnetization before one becomes dominant, ge- 
ometrical clusters do not compete. This picture 
is unfavorable for the use of geometrical clusters 
of Polyakov loops in gauge theories, for which the 
FK definition does not exist. 

The process of competitions between the 
largest FK clusters of different magnetization 
leads for the proper transition (h = 0) to a di- 
vergence of the equilibration time in the limit 
of infinite systems, an effect known in condensed 
matter physics 3 . Potts studies with an ex- 
ternal magnetic field show that a major slowing 
down effect survives for h 7^ 0. As the influence 
of an external magnetic field on the Potts model 
is similar to that of quarks on SU(3) gauge theory 
this indicates that the effect may be of relevance 
for QCD studies of the crossover region. 

For gauge theories a satisfactory cluster def- 
inition is presently not available. Nevertheless 
the underlying mechanism is expected to be the 
same as in the spin models. To study its influ- 
ence on the gluonic energy e and pressure p densi- 
ties, we like to calculate these quantities at times 

t — tmax • 

Let us first summarize the equilibrium proce- 
dure of Ref. |9I1 lj (in earlier work |12ll4j the 
pressure exhibited a non-physical behavior after 
the deconfining transition and the energy density 
approached the ideal gas limit too quickly due 
to the fact that anisotropy coefficients were cal- 
culated perturbatively). We denote expectation 
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values of spacelike plaquettes by P a and those 
involving one time link by P T . The energy den- 
sity and pressure can then be cast into the form 
(e+p)/T± = 



8N c N?g~ 2 



1 - yM a ) - CrO)] 



(P CT -P T )(11) 



and (e- 3p)/T 4 = 

\2N C N 4 [c a {a) - Cr(a)} [2P - (P CT + P T )] , (12) 

where Pq is the plaquette expectation value on a 
symmetric (T = 0) lattice and anisotropy coeffi- 
cients c a , T (a) are defined as follows: 



c<j, T (a) 



(13) 



They are related to the QCD /3-function in a sim- 
ple way 



dg~ 



da 



-2{c a {a) + c T {a)). 



(14) 



One needs to calculate the /3-function and 
anisotropy coefficients non-perturbatively to ob- 
tain meaningful results. In the range above the 
phase transition this is carried out in [3]. Using 
integral methods and Pade fits from JT] one finds 
the spatial anisotropy coefficient 



c a (a) = Co-(O) 



1 + dig 2 + d 2 g 4 
1 + d g 2 



(15) 



with d = -0.64907, d 1 = -0.61630 and d 2 = 
0.16965, where from perturbative calculations 
c CT (0) = 0.150702 iV c -0.146711/7V c . Thus, taking 
data for the /3-function from Table 3 of and 
evaluating c a (a) and c T (a) from Eqs. (|15|l and 
(|14l) we can calculate the energy density e and 
pressure p using Eqs. (|llfl and 112fl for arbitrary 
values of the coupling (3 above /3 C . 

To normalize to zero temperature, plaquette 
values from the symmetric N T — N a lattice are 
needed in Eq. I|12fl . For the quench this leaves 
us with the question whether we should perform 
the quench also on the symmetric lattice and 
use Po from the corresponding time evolution, or 
whether we should take the equilibrium value at 
(3 = 5.92 for Po- Fortunately, the empirical an- 
swer is that it does not matter. When one is far 
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Figure 4. SU(3) gluonic energy density: (a) 
with Po calculated from the time series after the 
quench and (b) using the equilibrium value of Pq. 

enough into the time evolution the results agree, 
see Fig. 0] for the time evolution on our N a = 16 
lattice. The reason is that one stays within the 
confined phase on the symmetric lattice. There- 
fore its equilibration after the quench is fast. 
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Figure 5. SU(3) gluonic energy density P(e) his- 
tograms: (a) with competing vacuum domains 
present and (b) after reaching equilibrium. 

Finally, we compare in Fig. [5] for the N a = 16 
lattice the gluonic energy distribution in equilib- 
rium at j3 = 5.92 with the one obtained after 
148 time steps. We find a shift towards lower 
gluonic energies and the width of the probabil- 
ity density is slightly broader for the time evolu- 
tion after the quench than in equilibrium. One 
also has to take into account that the geometry 
of relativistic heavy ion experiments experiments 
is reasonably approximated by N T /N,y — const, 
N a — ► oo, rather than by N T = const, N a — ► oo. 
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4. Summary and Conclusions 

Structure functions allow to identify the transi- 
tion scenario. For our quench from the disordered 
into the ordered phase of SU(3) lattice gauge the- 
ory we find spinodal decomposition. Relying on 
the linear theory of spinodal decomposition, we 
have calculated the critical mode k c . From it the 
Debye screening mass tod at temperature T is 
determined using phcnomenological arguments of 
Miller and Ogilvie [TU] . 

With increasing lattice size N a the time to 
reach the structure function maxima diverges. 
Relying on a study of Fortuin-Kasteleyn clusters 
in Potts models we assume that the reason 
is that vacuum domains of distinct Z3 trialities 
compete. These could be the relevant configura- 
tions after the heating quench in relativistic heavy 
ion experiments. We have initiated a study of the 
gluonic energy and pressure densities on such con- 
figurations. 

All our results rely on using a dissipative, non- 
relativistic time evolution, believed to be in the 
Glauber universality class. The hope is that 
the thus created non-equilibrium configurations 
may exhibit some features, which are in any dy- 
namics typical for the state of the system after 
the quench. This hope could get more credi- 
ble by studying a Minkowskian time evolution of 
Polyakov loops and finding similar features. Such 
a study appears to be possible \U>\ within a rel- 
ativistic Polyakov loop model which was intro- 
duced by Pisarski (IB"] . 
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